Enhancement of sound in chirped sonic crystals 
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We propose and experimentally demonstrate a novel mechanism of sound wave concentration 
based on soft reflections in chirped sonic crystals. The reported controlled field enhancement occurs 
at around particular (bright) planes in the crystal, and is related to a progressive slowing down of the 
sound wave as it propagates along the material. At these bright planes, a substantial concentration 
of the energy (with a local increase up to 20 times) was obtained for a linear chirp and for frequencies 
around the first band gap. A simple couple mode theory is proposed, that interprets and estimates 
the observed effects. The results are obtained for the case of sound waves and sonic crystals, however 
they are extendable to other type of waves in modulated host matter. 
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Manipulation and control of wave propagation, a prob- 
lem of fundamental interest, is at root of many applica- 
tions in different branches of science and technology. One 
important issue of wave manipulation is the localization 
and concentration (or local enhancement) of the wave 
energy. Artificial materials, and among them, artificial 
crystals are emerging as promising tools for manipulating 
wave propagation. In the case of sound waves considered 
here, such artificial periodic materials are called sonic 
crystals, structurally similar to photonic crystals in the 
field of optics. They are synthetic materials formed by 
a periodic distribution of elements or scatterers, whose 
properties (i.e., elasticity and density) differ from those of 
the host medium. This results in a periodic modulation 
of the acoustic properties of the medium at the scale of 
wavelength. The strong interest in these materials comes 
from their ability of manipulating the propagation of 
sound waves, due to their peculiar dispersive properties. 
A number of exotic and useful effects such as the forma- 
tion of band-gaps, 1,2 negative refraction^ birefraction^ 
self-collimation^ extraordinary transmission^ or local 
resonances^ have been so far demonstrated for sound 
waves. Utilizing these wave propagation effects, novel 
devices such as acoustic frequency filters ^ spatial (angu- 
lar) filters j2 lenses ^ or diodes^ have been proposed and 
demonstrated. 

We present here a novel wave propagation effect, con- 
sisting in specifically the wave energy concentration due 
to progressive decrease of the group velocity in chirped 
sonic crystals, in which the lattice constant, i.e. the dis- 
tance between scatterers in longitudinal (the wave prop- 
agation) direction, gradually changes along the propa- 
gation direction. We propose and demonstrate here a 
substantial increase of the wave intensity in controlled 
zones inside the crystal. Chirped (sometimes called 
graded) crystals have been introduced in optics 12 and 
acoustics^"— for different purposes, such as opening 
wide full band gaps or waveguiding of beams. An intrigu- 
ing phenomenon shown in chirped crystals is the smooth 
deflection of a light beam from the straight trajectory as 



it propagates through the crystal, the so-called mirage 
effect^. 

Another interesting effect reported recently is the so- 
called rainbow trapping effect, the dependence of the 
turning point position on the color of radiation. It has 
been predicted for one-dimensionally modulated chirped 
photonic structures^ and tapered optical and plasmonic 
waveguide s 18 i 19 . Rainbow trapping and wave enhance- 
ment are two different physical effects (the latter occurs 
even for monochromatic radiation), although they may 
occur simultaneously in chirped structures when the inci- 
dent radiation is broadband. In this letter, in addition to 
the extraordinary sound wave enhancement effect, which 
is the main result reported, we also present a sound rain- 
bow trapping effect for acoustic waves as a secondary 
result. 

Wave reflection from a band-gap in a chirped struc- 
ture is peculiar. The dispersion curves oj(k) at- and close 
to the band-edges, develop nearly horizontal segments, 
which corresponds to small or zero group velocity of the 
wave, since v g = duo/dh. The occurrence of the con- 
trolled sound enhancement requires that the crystal at 
the entrance plane be within the transparency range for 
the incoming wave, whose frequency is above the first 
band-gap, as shown in Fig. [TJ Note that here, not the 
wave frequency but the central (Bragg) frequency and 
the width of the band-gap is considered variable along 
the structure. Figures H{a) and DJb) show the variation 
of the band-gaps along the chirped sonic crystal, and the 
local dispersion curves at different depths, respectively. 
By local dispersion relation we mean the dispersion of an 
infinitely extended periodic crystal, for parameters (lat- 
tice constant, filling factor) corresponding to a particular 
depth of the chirped crystal. The wave entering into the 
crystal is gradually slowing down, as the local band-gaps 
are approaching the wave frequency in the course of prop- 
agation. Finally at a particular depth corresponding to 
the band-edge, the wave, literally speaking, stops, turns 
around, and starts propagating back. In other words it 
experiences a soft reflection. This effect is demonstrated 




FIG. 1. (Color online) (a) Dependence of the local band gaps 
on the local lattice constant along the chirped sonic crystal, 
(b) Band structure (local dispersion curves) evaluated at dif- 
ferent depths: at the entrance (top), at the middle (center) 
and at the exit (bottom) of the sonic crystal, (c) Intensity of 
the acoustic field calculated using Multiple Scattering tech- 
nique inside the chirped structure for the frequencies in (a) 
and (b). 

in Fig. Die), which shows the wave propagating through 
the crystal as obtained by numerical simulation using the 
multiple scattering theory approac h 20 ! 21 . The frequen- 
cies of the incident waves in the simulations correspond 
to local band-gaps at different depths. Figure QJc) evi- 
dences that the intensity of the wave increases substan- 
tially in the soft reflection area. Most importantly and 
in opposition to the case of perfectly periodic crystals 
(constant lattice period) in which only some discrete fre- 
quencies can be enhanced by the Fabry-Perot resonances, 
chirped crystals can localize the energy for a wide range 
of frequencies in a controlled way by the gradually change 
of the lattice constant. 

An experimental setup was designed to demonstrate 
the predicted extraordinary enhancement effect, and to 
obtain quantitative data of the acoustic field inside the 
structure. It consists in a two-dimensional sonic crystal 
with rectangular local symmetry, as illustrated in Fig. [2j 
made of acoustically rigid aluminum cylinders, of radius 
r = 2 cm, embedded in air. The spatial period is con- 
stant in transverse-to-propagation direction y, a y = 10 
cm, while a longitudinal chirp is introduced in the pe- 
riod along the propagation direction x. The adimensional 
chirp parameter is denned as a = (aj — aj+i ) / aj , where aj 
is the local longitudinal lattice constant at j-th layer. For 
our particular crystal case ao = 10 cm (initial period), 
ai3 = 4.8 cm (final period), and a gradient a = 0.055. 
The sign of the chirp can be either positive or negative, 
corresponding to lattice constant decreasing or increas- 
ing along the propagation direction. In case of identical 
scatterers, as used in our study, the filling fraction for 
the positive (negative) chirped structures increases (de- 
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FIG. 2. (Color online) Photographs of the experimental setup, 
(a) the grid of hanging points of the cylinders, (b) the chirped 
sonic crystal hunging vertically in the anechoic chamber. 



creases) in the propagation direction. This has a conse- 
quence of broadening of the local bandgap shown in Fig. 
Ha). 

The measurements were performed in an echo-free 
chamber with dimensions 8x6x3 m 3 , using the au- 
tomatized acquisition system 3DReAMS (3D Robotized 
e- Acoustic Measurement System) ^ This system enables 
the measurement of pressure fields along complicated tra- 
jectories as well as inside the crystals. Fig. [2fa) shows 
the grid of hanging points of the cylinders used to de- 
sign the chirped structure. Fig. [2fb) is a photograph 
showing the source (a loudspeaker) and the crystal in- 
side the chamber (for propagation directed downwards 
in Fig. [^a)). 

The experimental measurements are in excellent agree- 
ment with numerical calculations, as shown in Figs. [3] 
and [H We recorded the sound intensity profile along a 
fixed transverse position by translating the microphone 
along the x-axis through the void space between the rows 
of scatterers. In this way, we obtained two-dimensional 
space-frequency plots as shown in Fig. [3j from numerical 
(a) and experimental (b) data. White continuous lines 
mark the positions of the boundaries of the first band 
gap. Note the concentration of acoustic energy at posi- 
tions corresponding to just before the upper band edge. 

Figure [3Ja) shows the axial distributions obtained ex- 
perimentally (dots) and theoretically (continuous lines) 
for three particular frequencies. In both cases small-scale 
fringes are observed, corresponding to the local Bloch 
mode, as well as a large-scale oscillations or envelope 
(dashed line) of the Bloch mode, to be discussed below. 
Figure [U(b) represents the theoretical calculation of the 
position of the maximum value of concentration of acous- 
tic energy inside the crystal depending on the frequency. 
In correspondence with the results in Fig. [TJ the position 
of the maximal energy concentration shifts deeper into 
the bulk of the structure, as the frequency is increased. 
Note also that, since the incident amplitude was normal- 
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FIG. 3. (Color online) Sound distribution inside the crystal, 
along the x-axis, for varying frequencies (a) Multiple Scatter- 
ing simulation and (b) experiment. White continuous lines 
show the spatially varying edges of local band-gaps. Color 
scale represents the values of the acoustic intensity. 



ized to unity, at the maximum value, the intensity has 
been recorded up to around 20 times higher than inci- 
dent. For usual reflection between two different homoge- 
neous media or from a purely band-gap material in the 
range of the band-gap, only an increase of 4 times of the 
local intensity is possible (as the interference pattern is 
formed from forward and fully reflected backward wave) . 
For the case of periodic structures the wave penetrates 
into the reflecting material evanescently, i.e. with expo- 
nential decay22 and never shows an increase of intensity. 

The increase of the intensity field observed in Fig. 2J 
can be understood from a coupling between the forward 
and backward waves inside the chirped crystal. In or- 
der to interpret the results, we propose a simple coupled 
mode analytical theory for the propagation of acoustic 
plane waves inside a one-dimensional chirped crystal (in 
optics also known as chirped mirror). This dimension 
reduction is possible because the first band gap in TX 
direction essentially appears due to a resonant coupling 
between the forward and the backward waves. The con- 
tributions of the wave components propagating to trans- 
verse directions are negligible. Assuming that the full 
pressure field consists of forward and backward propa- 
gating waves, P = A(x)e zkx - lujt + B(x)e~ zkx - ZUJt + c.c. 
the following coupled amplitude equations can be sys- 
tematically obtained from wave equations, 



dA 

dx 



dB 
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a(x) 



a(x) 
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where s is the back-reflection coefficient by one row 
of scatterers, a(x) the variable longitudinal period and 
Aq(x) = 2tt/X — ir/a(x) is the detuning from the Bragg 
frequency. 

From the numerical study of the scattering by only 
one row of the structure, we estimate that the back re- 
flected intensity is around 40% of the incident, so s c± 0.4. 
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FIG. 4. (Color online) (a) Numerical simulation results (con- 
tinuous line) and experimental results (dots) for the acoustic 
intensity at the central section inside the crystal for the three 
frequencies: 2500 Hz, 2600 Hz, and 2700 Hz. The shaded 
(yellow) rectangle denotes the area covered by the crystal, 
(b) Position of the maximum value of concentration of energy 
inside the crystal depending on the frequency of the incident 



The same numerical study reveals that the scattering 
into field components propagating at transverse direc- 
tion is only around 5%, which justifies the followed one- 
dimensional approach, neglecting transverse modulations 
in the vicinity of the first bandgap. We notice that 
the detuning from the Bragg resonance Aq(x) is a func- 
tion of the longitudinal position x for chirped crystals. 
Recall that in our study the chirp is linear, given by 
a(x) = do + a(x — xq). 

Equations (pQ) can be rewritten in canonical form as 



d 2 A 
dX 2 



ie(x) 



dA 
dX 



A, 



(2) 



where the space scaling dX = dx^J~s/a(x) was chosen 
to make the normalized coupling coefficient unity, and 
e(x) = 2d(XAq(X)/dX is the normalized detuning from 
the Bragg frequency. 

The wave, roughly speaking, reflects from the bandgap, 
i.e. from the position Xo corresponding to the Bragg fre- 
quency, with e(Xrj) = 0. In general (for arbitrary chirp) 
Eq. ([2]) cannot be solved analytically. However, in a 
simple case when the normalized detuning varies linearly 
around zero e(X) = ei(X—Xo), Eq. (j2]) has an analytical 
solution in the form 



A(X)= Cl H t/ei (Xyfe/2), 



(3) 



where H n is the Her mite polynomial of imaginary order. 
The counter-propagating field obeys a sim ilar ex pression. 
The integration constant c\ = H % / tl {Xp y it\ /2) is deter- 
mined by the boundary conditions, by imposing that the 
amplitude of the forward wave at the front face X = Xp 
equals unity. e\ = de(X)/dX\x=x F or > m terms of initial 
variables, e\ = Ana/ which estimated for experimental 
parameters results e± = 3. 




FIG. 5. (Color online) Solutions of (3), i.e. the intensity of 
the forward (continuous line) and backward (dashed line) field 
along the chirped structure, (a) (with ) corresponds to exper- 
imental configuration, (b) (with ) is shown for comparison, to 
illustrate soft reflections for a substantially smaller chirp. The 
vertical dashed lines indicate the center of the bandgap, and 
the shaded areas the bandgap itself. 



In Fig. [5]we present the amplitude of the acoustic in- 
tensity of the forward and backward waves for linearly 
chirped crystals as follows from Eq. (|3j). The acoustic 
field is nearly exponential in the bandgap, and oscillatory 
in front of it. The oscillations, with the period and am- 
plitude increasing as the wave approaches the band-gap, 
are large-scale oscillations, which originate from the en- 
ergy exchange between the forward and backward waves. 
These large-scale oscillations correspond to oscillations of 
the envelope of the Bloch modes observed in Fig. 2J and 
are not due to conditions imposed at the entrance of the 
sonic crystal, e.g. some possible impedance mismatch. 

The controlled field enhancement effect is clearly vis- 
ible again in Fig. [5] From the analytical estimations in 
Eqs. (P)-© and from Fig. [5] it follows that for maximal 
field enhancement of the wave intensity, the chirp must 
be as small as possible. For e\ = 0.3 the maximum field 
enhancement could be around 6 times (in terms of inten- 
sities) if one compares the maximal and minimal values 
of the plot in Fig. [5(b). In order to realize such enhance- 
ment the entrance to the sonic crystal must be placed to 
correspond to the deepest minimum of the solution (j3j), 
in this particular case at around the point X =~ 10. 
For the parameters of the experiment, e\ = 0.3, the en- 
hancement of more than two times is predicted in this 
simplified approach. Also, as Fig. [5(a) shows, a small 
portion of radiation is transmitted, i.e. leaks through the 
band gap. Such tunneling, analogous to Landau-Zenner 
tunneling, is due to slightly too fast (no more adiabatic) 
chirp. 

Concluding, in this Letter we have predicted and ex- 
perimentally demonstrated a novel mechanism for ex- 
traordinary sound field enhancement in a chirped crys- 
tals, specifically in chirped sonic crystals. The acoustic 
wave energy can be selectively concentrated at particular 
depth of the crystal depending on the frequency and on 



the parameters of the structure. At these bright planes, 
a substantial increase of the energy was recorded for lin- 
ear chirp and for frequencies around the first gap along 
the X direction of structure. The experimental study 
was performed in a macroscopic sonic crystal irradiated 
by acoustic waves in audible regime, where the measure- 
ments in the interior of the crystal are possible. 

We note, that the reported enhancement of the field is 
of different nature from the enhancement of the fields in 
resonant systems, e.g. in Fabry-Perot resonators, or in 
finite size periodic (unchirped) structures. The field in- 
side the Fabry-Perot-like systems increases by resonator 
quality factor Q, however the enhancement occurs for a 
narrow frequency range (proportional to Q _1 ). Therefore 
in average the intensity of broad frequency band radia- 
tion does not increase at all in Fabry-Perot-like resonant 
systems. The chirped sonic crystal studied here is not 
a resonant system, and the field enhancement occurs for 
broad frequency range (just the position of field localiza- 
tion is frequency dependent). 

A simple couple mode theory is proposed, that inter- 
prets and estimates the observed effects. The theory is 
based on generic equations describing wave propagation 
in periodically modulated media. The presented results 
have a universal character, i.e. valid irrespective to the 
nature of the waves. This implies that the reported ef- 
fect, the controlled enhancement of the wave energy can 
be observed e.g. for optical waves, microwaves, or per- 
haps some exotic waves as surface plasmon waves. 

In the field of acoustics the results are independent of 
the spatial scale of the structure, and in principle the phe- 
nomenon could be scaled-down and observed in micro- 
or nano-scale phononic (so called hypersonic) crystals^. 
At these scales sound waves are described in terms of 
phonons, and the ideas presented in this work could 
find application for heat management in acoustical or 
acousto-optical devices. Recent works in this direction 
show indeed that manipulation of phonon dispersion 
properties can allow thermal transport control^!. Gen- 
erally, the effect of wave energy concentration demon- 
strated in present work, opens a possibility of increasing 
the efficiency of detectors and absorbers, both in acous- 
tics and optics, since slow phonons and photons can be 
absorbed and harvested with a higher probability. 



ACKNOWLEDGMENTS 

The work was supported by Spanish Ministry of 
Science and Innovation and European Union FEDER 
through projects FIS2011-29731-C02-01 and -02, also 
MAT2009-09438. V.R.G. is grateful for the support of 
post-doctoral contracts of the UPV CEI-01-11. K.S. ac- 
knowledges the grant of UPV PAID-02-01. We acknowl- 
edge the Centro de Tecnologias Fisicas: Acustica, Materi- 
ales y Astrofisica and the Sonic Crystal Technologies Re- 
search Group of the Universitat Politecnica de Valencia 
for the use of the anechoic chamber and the 3DReAMS 



5 



respectively. 



1 M. Kushwaha, P. Halevi, L. Dobrzynski, and B. Djafari- 
Rouhani, Phys. Rev. Lett. 71, 2022 (1993). 

2 R. Martmez-Sala, J. Sancho, J. V. Sanchez, V. Gomez, 
J. Llinares, and F. Meseguer, Nature 378, 241 (1995). 

3 X. Zhang and Z. Liu, Appl. Phys. Lett 85, 341 (2004). 

4 M.-H. Lu, C. Zhang, L. Feng, J. Zhao, Y.-F. Chen, Y.-W. 
Mao, Y.-Y. Zhu, S.-N. Zhu, and N.-B. Ming, Nature Mat. 
6, 744 (2007). 

5 V. Espinosa, V. J. Sanchez-Morcillo, K. Staliunas, I. Perez- 
Arjona, and J. Redondo, Phys. Rev. B 76, 140302(R) 
(2007). 

6 Y. Zhou, M.-H. Lu, L. Feng, X. Ni, Y.-F. Chen, Y.-Y. Zhu, 
S.-N. Zhu, and N.-B. Ming, Phys. Rev. Lett. 104, 164301 
(2010). 

7 Z. Liu, X. Zhang, Y. Mao, Y. Zhu, Z. Yang, C. Chan, and 
P.Sheng, Science 289, 1734 (2000). 

8 A. Khelif, P. A. Deymier, B. Djafari-Rouhani, J. O. 
Vasseur, and L. Dobrzynski, J. Appl. Phys., Vol. 94, 1308 
(2003). 

9 R. Pico, V. Sanchez-Morcillo, I. Perez- Arjona, and 
K. Staliunas, Appl. Acoust. 73, 302 (2012). 

10 F. Cervera, L. Sanchis, J. Snchez-Prez, R. Martnez-Sala, 
C. Rubio, and F. Meseguer, Phys. Rev Lett. 88, 023902 
(2002). 

11 X. Li, X. Ni, L. Feng, M. Lu, C. He, and Y. Chen, Phys. 



Rev. Lett. 106, 084301 (2011). 

12 E. Cassan, C. K.-V. D. Caer, D. Marris-Morini, and 
L. Vivien, J. Lightwave Tech. 29, 1937 (2011). 

13 M. Kushwaha, B. Djafari-Rouhani, L-Dobrynski, and 
J. Vasseur, Eur. Phys. J. B 3, 155 (1998). 

14 I. E. Psarobas and M. M. Sigalas, Phys. Rev. B 66, 052302 
(2002). 

15 L. Wu and L. Chen, J. Appl. Phys. 110, 114507 (2011). 

16 E. Centeno, D. Cassagne, and J. -P. Albert, Phys. Rev. B 
73, 235119 (2006). 

17 Y. Shen, J. Fu, and G. Yu, Phys. Lett. A 375, 3801 (2011). 

18 V. N. Smolyaninova, I. I. Smolyaninov, A. V. Kildishev, 
and V. M. Shalaev, Appl. Phys. Lett. 96, 211121 (2010). 

19 M. Jang and H. Atwater, Phys. Rev. Lett. 107, 207401 
(2011). 

20 P. Martin, Multiple Scattering. Interaction of Time- 
Harmonic Waves with N Obstcles (Cambirdge University 
Press, UK, 2006). 

21 Y. Y. Chen and Z. Ye, Phys. Rev. E 64, 036616 (2001). 

22 V. Romero- Garcia, J. Sanchez-Perez, and L. Garcia-Raffi, 
J. Appl. Phys. 108, 044907 (2010). 

23 T. Gorishnyy, C. K. Ullal, M. Maldovan, G. Fytas, and 
E. L. Thomas, Phys. Rev. Lett. 94, 115501 (2005). 

24 P. Hopkins, C. Reinke, M. Su, R. O. Ill, E. Shaner, Z. Lese- 
man, J. Serrano, L. Phinney, and I. El-Kady, Nano Lett. 
11, 107 (2011). 



